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1. INTRODUCTION 

Neutron stars are the densest stable stellar objects in the present day Universe. 
For the typical masses M ~ 1 — 2 Mq, they are expected to have radii R ~ 
10 km, so that their mean density p ^ 10^^ g cm~^, significantly higher than 
the normal nuclear density, po — 2.5 x lO^'^g cm~'^, characteristic of interiors 
of heavy atomic nuclei. Neutron stars are formed as very hot objects, as an 
outcome of gravitational collapse of massive stellar cores at the endpoint of their 
thermonuclear evolution, or - under specific conditions - in the gravitational 
collapse of mass accreting white dwarfs. Newly born neutron stars are very 
hot, with Tinitiai ^ 10" K (which corresponds to A:bT - 10 MeV). They cool 
rapidly due to neutrino emission from their interiors, so that after one year their 
internal temperature falls below 10^ K. At such temperatures, matter in the 
neutron star interior is strongly degenerate (due to its huge density), so that 
the effects of temperature on the structure of neutron star, older than one year, 
can be neglected. Neutron stars are observed as radio pulsars. X-ray pulsars, 
and as X-ray bursters. Since they contain matter of the density ranging from 
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a few g cm~'^ in the atmosphere to >10^^ g cni~^ at the center, they should 
be considered as 'cosmic laboratories' for superdense matter. Observations of 
neutron stars yield constraints on the theory of dense matter which was used for 
constructing neutron star models (an introduction to physics and astrophysics 
of neutron stars can be found in 

Neutron stars are massive (M 1 — 2 Mq), compact {2GM/Rc^ ^ 0.3), and 
rapidly rotating (period of rotation P ~ 1 ms for millisecond pulsars) - and 
therefore they can be considered as potentially promising sources of continuous 
gravitational radiation. The shortest observed period of rotation (for radio 
pulsar PSR 1937+21) P = 1.56 ms, which corresponds, for R ~ 10 km, to 
the equatorial velocity i^equator ~ 0.1c. The expected typical moment of inertia 
of neutron star is / ~ 10^^ g cm^. Crucial for the emission of continuous 
gravitational radiation is the time dependence (in the laboratory frame) of the 
quadrupole moment tensor, Qij, of the stellar mass distribution. Let us denote 
the amplitude of the time dependent part of Qij by /eoff , where CeS describes 
the deviations from the axial symmetry (with respect to the rotation axis), 
and is ususally referred to as effective tri-axiality of neutron star [Q. The mean 
power of the gravitational radiation, emitted by rotating neutron star with 
non- vanishing Ceff, is proportional to e^g, and strongly increases with rotation 
frequency, Egr oc /2e2gP-6 (|. 

The non- vanishing value of eoff appears only under specific conditions. It can 
be sustained due to a partially solid nature of neutron star (a completely fluid, 
stably rotating body is symmetric with respect to rotation axis). As discussed 
in Sections 6-8, the existence of the solid interiors of neutron stars can lead to 
a non-zero e^s in two cases: first, due to the appearance of genuine tri-axiality, 
resulting from the presence of 'mountains' on the neutron star crust, or due to 
the neutron star precession. Two other possibilities, not related to the solid 
nature of neutron stars, and therefore not discussed in these lectures, should 
be mentioned. The first possibility corresponds to strong magnetic field, with 
symmetry axis different from the rotation axis. However, such magnetic field 
should be extremely strong (i3>10^^ G), and should have quite special spatial 
structure to maximize eeff (see the lecture of S. Bonazzola and E. Gourgoulhon 
in this volume). The second possibility, valid also for a completely fluid star, is 
connected with appearance of non-axisymmetric, gravitational radiation reac- 
tion driven instabilities, |^ (see also the lecture of R. Wagoner in this volume), 
or with the shape instabilities, discussed in this volume in the lecture of S. 
Bonazzola and E. Gourgoulhon. However, both gravitational radiation reac- 
tion and shape instabilities are expected for a rotation frequency f2 = 2tt / P 
close to the mass-shedding limit, f2 ~ J^mSi while observed 17 's are significantly 
lower than JIms- 

An additional attractive feature of a (partially) solid neutron star with non- 
zero eoff is connected with the remarquable stability of pulsar rotation. This 
could allow for long integration times (years) of periodic gravitational radiation 
signal from a pulsar with non-zero Ceff. 

The existence of a solid crust of neutron star is well established on the theo- 
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retical ground. Theory predicts the sohd character of dense, cool (degenerate) 
matter for the density below ph — lO^"* g cm~^. The existence of the solid 
crust is confirmed in a rather convincing way by the observations of glitches in 
the pulsar timing 0|. On the contrary, the possibility of existence of the solid 
cores of neutron stars, with density 10^^ g cm~'^, is based on rather exotic 
models of superdense matter, and should be considered as interesting but very 
uncertain. 

2. FORMATION AND COMPOSITION OF THE SOLID CRUST 

In the present section we will restrict ourselves to the outer layer of the neutron 
star crust, with density up to a few times 10^'^ g cm~^. In this density regime, 
the volume occupied by the atomic nuclei constitutes a small fraction of the 
total volume and atomic nuclei, present in the matter, can be treated as spher- 
ical. The properties of the bottom layers of the crust, with p ~ 10"'^'* g cm"'', 
within which the assumption of sphericity of nuclei may break down, will be 
discussed in Section 3. 

2.1. Formation of the crust: high temperature scenario 

The initial temperature of the interior of neutron star formed in gravitational 
collapse is Tjnitiai ^ 10^^ K, which corresponds to the typical energy of thermal 
motion of matter constituents ^ 10 MeV. At such temperature matter is a 
multicomponent plasma of atomic nuclei (A, Z), alpha particles a, neutrons n, 
protons p, and electrons and positrons e~ , e^. The baryon density of matter 
(number of baryons - in our cases nucleons - in a unit volume) is given by 
nh — n-n + np + + tiaA^ where rii is the number density of a species 'i', and 
where we assumed, for simplicity, that only one heavy nuclid (j4, Z) is present 
in the matter (this is a rather good approximation, see |^). Charge neutrality 
of matter requires that + + 2na + Zua = J^e- • Under the conditions 
prevailing in the interior of a newly born neutron star, matter can be assumed 
to be in complete nuclear equilibrium (the rate of all reactions is much higher 
than the rate of evolution of neutron star interior). At given Tib and T, the 
composition of matter is obtained from the minimization of the free energy 
density / = e — Ts, where e is the internal energy density, and s is entropy 
per unit volume. Neutron star cools rapidly, due to neutrino losses, and after 
one year internal temperature is less than 10® K. A standard assumption is, 
that during the cooling process matter was always close to the state of nuclear 
equilibrium, with composition always given by the minimum of /. However, at 
T < 10® K thermal contribution to / is negligible (for p > 10^ g cm"''), and the 
T — approximation (ground state of matter) is excellent. The composition of 
matter can be then calculated from the requirement e = minimum, at fixed ?ib, 
and under the condition of charge neutrality. Such a procedure of determining 
the ground state of dense matter enables one to calculate the composition and 
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equation of state of the outer envelope of neutron star, formed from the initial, 
very hot state. The ground state of matter turns out to be a body-centered- 
cubic (bcc) lattice of nuclei {A,Z), immersed in the uniform background of 
electrons. The neutron fraction in nuclei, N/Z, increases with density, and 
above neutron drip density, pnd = 4 — 6 x 10^^ g cm~'^, some of neutrons form 
a gas outside nuclei, so that for p > pnd ground state of matter is a bcc crystal 
of nuclei, embedded in a neutron gas, and permeated by a uniform background 
of electron gas. 

2.2. Formation of the crust: low temperature scenario 

Many neutron stars are (or were in the past) members of close binary systems, 
in which the second member is (was) an evolved, large radius (and therefore low 
surface gravity) star. Due to its huge gravitational field, neutron star accretes 
matter from its companion. Accretion of matter increases the mass of neutron 
star, leads to formation of a new outer mantle, and influences the tempera- 
ture distribution within the neutron star interior. After the quasi-stationary 
thermal state of the neutron star interior has been reached, the internal tem- 
perature profile can be obtained from the balance of energy gains and losses. 
The energy gains are due to heating by infalling matter and to work done by 
gravitational forces compressing the stellar interior. The heat loss is due to 
photon emission from the neutron star photosphere and neutrino losses from 
the neutron star interior. The calculations, based on realistic neutron star mod- 
els, lead to conclusion, that for typical accretion rates in close binary systems, 
10~^^ Mq/y < M < 10~^ Mq/y, internal temperature is always lower than 
10^ K 1^, ^, Ij. At such temperatures the assumption of full thermodynamic 
equilibrium during the formation of the new neutron star crust is not valid: 
the energy of thermal motion of nuclei, -Ethcrm = f fcer, is so low compared to 
typical Coulomb barriers, -Ecoui, that nuclear reactions involving charged par- 
ticles are practically suppressed (blocked) (on the timescale of the neutron star 
evolution). The pressure and density in a given shell of the crust of accreting 
neutron star increase, as the accretion proceeds. The local value of electron 
chemical potential (Fermi energy) increases according to /ig oc (pZ/A)^/"^, for 
p > 10^ g cm^"^. However, the composition of the matter can change only via 
those reactions between matter constituents, which are energetically possible 
under prevailing low temperature conditions. The only reactions, which can 
proceed under conditions prevailing in the neutron star crust are thus: electron 
capture on nuclei, 

(A, Z) e- — > (A, Z - 1) + i^e , (1) 
neutron emission (if triggered by electron captures) , 

{A,Z)+e- — > {A-k,Z -l) + iye + k-n , (2) 

as well as neutron absorption, and pycnonuclear fusion of neighbouring nuclei, 
resulting from quantum tunnelling of the Coulomb barrier due to quantum 
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zero-point motion around the crystal lattice sites, 

iA,Z) + {A,Z) — y {2A-k,Z) + k-n . (3) 

The pycnonuclear regime of fusion (Greek pycnos - means dense) is to be con- 
trasted with thermonuclear one; the kinetic energy of nuclei results there from 
the quantum effect of localization, which at sufficiently high density dominates 
over the thermal one: Ezcm point ^ -^thermal- The value of Ezcm point increases 
with increasing density, which results in a strong density dependence of the 
rate of pycnonuclear fusion jl] . 

The crust of an accreting neutron star has a characteristic, time dependent 
layer structure. Let us assume, that during time t the accretion of matter onto 
the neutron star surface took place at constant rate M. After time t, the outer 
layer of neutron star of mass AM = Maccr = Mt has thus been replaced by the 
accreted matter. This new shell of mass AM compresses the 'old' crust, laying 
below it. Each layer of the 'old crust' undergoes compression, characterized by 
the compression factor, equal to the ratio of its actual density, p (at time i), to 
its initial, original density, pi (before accretion). Let us characterize the original 
(initial) position of the crust layer by the mass of the layers laying above it, 
AMi, and let us denote the corresponding compression factor by C(AMi) = 
p/pi- Notice, that the pressure at the bottom of the layer is proportional to 
the mass laying above it. Consider the case of Maccr = 5 x 10"'' Mq, which 
corresponds to t = 5 x 10^ y at accretion rate M — 10^^^° Mq/j. In what 
follows, we will quote values obtained for a specific model of neutron star, of 
mass 1.4 Mq, with liquid core (p > ph) calculated for the BJI equation of state 
For a layer, laying originally at the neutron star surface (for this layer 
AMi = 0), one gets then C = lO'^'^: the density of this layer after accretion 
of AMaccr is ~ 10^^ g cm~'^. However, due to huge pressure gradient in the 
surface layers, the value of C decreases rapidly with increasing AM;. We get 
C(AMi = 5 X 10-* Mq) ~ 2 (i.e., the density of this layer increase from 
10" g cm-3 to 2 X 10" g cm-3). Finally, we get C(AMi = 10"^ Mq) ~ 1. 
The compression of the deeper layers of neutron star crust can be neglected, 
and so their composition coincides with the original one. 

Let us summarize our results. If neutron star accreted Maccr, then the compo- 
sition of its outer layer of mass Maccr may be expected to deviate widely from 
the original (ground state) one: this layer of neutron star crust was formed 
under conditions differing strongly from the 'high temperature' scenario, de- 
scribed in the preceding section. Its actual composition can be obtained by 
following the evolution of matter, which was originally a hydrogen rich plasma 
lost by the neutron star companion. Such a program was fulfilled in 0, ||. The 
differences with respect to the ground state decrease gradually with increasing 
density, and for the crust layers laying below the outer envelope of the mass 
2Maccr the approximation of the ground state composition is very good. 
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2.3. Composition of the crust, melting temperature 

Composition of the crust is important for determining its pliysical properties, 
such as melting temperature, elastic parameters, and thermal and electrical 
conductivity. Within standard approximation, the composition of neutron star 
crust will be determined by giving the values of A and Z of the nuclid present in 
the matter and, above the neutron drip point, also the mass fraction contained 
in the neutron gas outside nuclei, X„. These parameters, for the ground state 
of the crust, and for the case of accreted crust, are shown in Figs. 1, 2. As we 
see from Fig. 1, the values of Z and A in accreted crust are significantly smaller, 
than those characteristic of the ground state of matter. The differences increase 
with increasing density, and for p > 10^^ g cm~'^, Z and A in the accreted crust 
are 2-3 times smaller than the ground state one. As seen from Fig. 2, neutron 
drip in an accreted crust occurs at a somewhat higher density, than in the case 
of the ground state of dense matter, but the values of Xn at higher densities 
are not dramatically different. In both cases, matter undergoes neutronization: 
the proton mass fraction decreases as the density increases. 

One notices characteristic features of the density dependence of A and Z, 
represented in Fig. 1. They can be understood in terms of the general proper- 
ties of nuclei. Only even-even nuclei (even numbers of neutrons and protons) 
are present - which results from the action of the pairing forces between nu- 
cleons. In the case of the ground state of dense matter, strong shell effects 
are clearly seen: this results from the fact, that closed-shells nuclei are par- 
ticularly strongly bound. For density below 10^ g cm""^, increasingly neutron 
rich Ni isotopes are present, with Z = 28 (closed proton shell). Then, up to 
the neutron drip point, nuclids with closed neutron shells iV = 50 and = 82 
are present in the ground state of matter. Finally, above neutron drip, closed 
Z = AO proton subshell and Z = 50 proton shell are persistent, while neutron 
shell eff'ects are no longer visible (which is due to the presence of the outer, un- 
bound neutron gas). Let us notice, that up to the density p ~ 10^^ g cm~^ the 
composition of the ground state of matter is determined by the experimental 
data 1^ ! In the case of accreted matter, one starts with ^^Fe shell (produced 
from the electron captures on ^^Ni, obtained at the considered accretion rate 
by the explosive burning of helium formed by steady burning of accreted, orig- 
inally hydrogen rich plasma). Then one has to follow the evolution of this shell 
during low temperature compression up to ~ 10^^ g cm^'^ 

Up to now, we neglected the effect of finite temperature of neutron star 
crust. While at T < 10^ K the effect of T on the composition of matter is 
negligible, the value of T decides whether, at a given density, dense plasma 
is in the liquid or solid (crystal) phase. Crystallization of the plasma is due 
to Coulomb forces between its constituents, which lead to localization of ions 
(nuclei). The importance of Coulomb energy, as compared to thermal energy 
of nuclei, is measured by the Coulomb coupling parameter. 
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Fig. 1. — The values of Z (upper panel) and A (lower panel) for the nuclids present 
in dense matter for the ground state (solid line: based on [5] below neutron drip point 
and [4] above neutron drip) and for the accreted crust (dashed line: based on [6]). 

where a is the radius of the sphere of the volume 'per one nucleus', l/n^ 
(volume of the Wigner-Seitz (W-S) ceU, yccii IT])- The many-body calculations 
of the thermodynamic properties of dense one-component plasma show that, 
at fixed composition, solid phase of the plasma (bcc crystal) is energetically 
preferred over the liquid one for F > ~ 178 ||ll|. This corresponds to the 
melting temperature 

\_ 

= 1.32 X 10^ (^] ' K , (5) 

where Aceii is the number of nucleons in the Wigner-Seitz cell, and pg — 

(^) The Wigner-Seitz cell is a useful object in the plasma physics; in our ap- 
proximation it is a sphere containing one nucleus and Z electrons, to be electrically 
neutral. To a very good approximation, the energy of a plasma is a sum of energies 
of Wigner-Seitz cells. 
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Fig. 2. — Mass fraction contained in the neutron gas outside nuclei, X„, for the 
ground state of dense matter (solid line) and accreted crust (dotted line) 



p/(10^ g cm ■^). For p < pnd, we have Accu = A, while for p > pnd the 
total number of nucleons in the W-S cell is ^coii = A + UnVccn- 

The melting curve T = T^{p) separates the liquid and solid regions in the 
density - temperature plane. The melting curves for the ground state of dense 
matter, and for the accreted crust, arc shown in Fig. 3. The difference between 
the two melting curves are significant for p > 10^ g cm""^, and above neutron 
drip point T^{ground state) » T^{accreted crust). This reflects strong Z- 
depcndence of melting temperature, combined with rather weak dependence 
on Acei\ and p. 



3. THE BOTTOM LAYERS OF THE CRUST: EXOTIC NUCLEAR 
SHAPES AND TOPOLOGIES 

For p > pnd nuclei can be considered as droplets of 'nuclear matter', of nearly 
constant density pdrop — Po = 2.5 x 10^^ g cm~'^, immersed in a less dense 
neutron gas. So, above neutron drip density, nucleons are present in two coex- 
isting phases: gaseous (neutron gas, or G-phase) and liquid (nuclear matter, or 
L-phase) . As long as the fraction of the volume occupied by nuclei is a small 
fraction of the total volume of the system, droplets of nuclear matter can be 
treated as spherical - such a shape minimizes the surface energy of the droplets. 
At T = 0, the energy density at a given baryon density rib can be decomposed 
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Fig. 3. — Melting temperature versus density, for the ground state composition of 
dense matter (solid line) and for the accreted crust (dashed line) 



as 

e(nb, parameters) = evoi.A + Cgurf + ecoui + eg + e„ , (6) 

where emi.A is the bulk (volume) part of the contribution from nuclei (droplets) , 
Csurf is the contribution, resulting from the existence of the interface between 
the neutron gas and nuclear matter (surface energy), and ecoui is the contribu- 
tion from Coulomb interactions of charged constituents of matter. Finally, Ce 
is the energy density of electron gas, which permeates both phases of nucleonic 
matter, and e„ is the contribution of neutron gas outside nuclei (droplets). 
The description of nuclei embedded in neutron gas is usually done using com- 
pressible liquid drop model of atomic nuclei. Within this model, the problem 
of determination of the ground state of matter at some fixed n^, is reduced to 
minimization of e(nb, parameters) with respect to all parameters, under the 
condition of charge neutrality. Let us notice, that the parameters of 'droplets' 
are modified with respect to the zero pressure case (characteristic of terrestrial 
nuclei) due to pressure exerted by neutron gas (which results in compression 
of nuclei); the presence of neutron gas influences also the surface energy term 
in Eq. §). 
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The fraction of volume occupied by the denser L-phase increases with matter 
density: however, even at p = 10^'^ g cm~^ we have Vl/F = 2 x 10"^ Q, 
and the approximation of spherical nuclei can be expected to be valid. The 
situation changes, however, when we consider matter of the density close to po, 
where Vl constitutes a sizable fraction of V: then, the spherical droplets may 
no longer minimize the energy density of the system. Such a possibility was 
considered long time ago by Baym, Bethe and Pethick who showed that 
at Vi^/V = 1/2 nuclei turn 'inside out', and at higher Vi,/V spherical bubbles 
of neutron gas in nuclear matter are energetically preferred. 

The problem of actual nuclear shapes for p^/S < p < po was recently recon- 
sidered by Lorentz et al. fl^ . Using the compressible liquid drop model, these 
authors optimized the value of e(nb, parameters, shape) with respect to the 
shape of the interface between the L and G-phase. The shapes they considered 
were: in three dimensions (3-D) - spherical droplets of L-phase in homoge- 
neous G-phase, and spherical bubbles of G-phase in homogeneous L-phase; in 
2-D - cylindrical rods of L-phase in homogeneous G-phase, and cylindrical tubes 
filled with G-phase in homogeneous L-phase. Finally, considered 1-D shapes 
were alternating plates of G and L-phases. The sequence of (first-order) phase 



transitions, obtained in |13 for a specific compressible liquid drop model of 
nuclei, starts a,t p ~ 1.06 x 10^"^ g cm~^, where Vt = 0.09, X„ = 0.76 and 
Xp = 0.03, by a transition to 'rods'. Then, with increasing density, transitions 
from 'rods' to 'plates', from 'plates' to 'tubes', and from 'tubes' to 'bubbles' 
take place. Finally, at ph = 1-6 x lO^'' g cm~^ matter becomes a homogeneous 
mixture of neutrons, protons and electrons (Vl = V, Xp = 0.03, X„ = 0.97). 
The value of ph corresponds to the bottom of the solid crust as determined in 



13 . This value depends somewhat on the assumed parameters of the com- 
pressible liquid drop model, used for the description of the the L-phase, and 
on the n — n interaction model, relevant for the description of the G-phase. 



4. CRUST CONTRIBUTIONS TO STELLAR MASS AND MO- 
MENT OF INERTIA 

In order to calculate the contribution of neutron star crust to neutron star 
parameters, such as mass and moment of inertia, let us consider the equation 
of hydrostatic equilibrium for a non-rotating neutron star, 

dp Grnpf^^4nrfp\ 



dr r 

dT = '""^'^ ^=l^"7^j ■ 

For sufficiently massive neutron star (A/>1.2 M©), the crust is thin, and con- 
tains only a small fraction of the mass of neutron star. For a shell at radius 
r within the crust, both {R — r)/R and (M — m)/M are typically few times 
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10"^. Moreover, pressure within the crust is smaU compared to energy density, 
vl p(? < 10"^. AU these conditions enable us to rewrite Eq.(Q) within the crust 
in the approximate form, 

dp GMpKjr) dm ^ 

T ^ 2 ' — ^Anr p , 8 

ar ar 

which can be rewritten, with our approximation, as 

dp GMpk{R) 



dm 47ri?2 



(9) 



Integration of the above equation imphes, that the mass of the layer of the 
crust is proportional to the pressure at its bottom, 

" Gum ■ 

Similar approximations, applied to the moment of inertia for slow, rigid rota- 
tion, lead to the formula [|l3j 

A/(pbot) ^ -AM(pbot)i?' f 1 - A(i?) , (11) 

where / is the total moment of inertia of the star. 

Combining these equations with equation of state for the crust, one can see 
that the innermost layer of the crust, containing 'exotic nuclei', discussed in 
the preceding section, yields as much as half of the total crust contributions to 
M and / ! For a 1.4 Mq neutron star, /crust and Afcrust constitute ~ 1 % of / 
and Af , respectively (this agrees well with exact results reported in [^). Both 
A/crust/AT and /crust// decrease with increasing M. 

The value of ph, quoted in this lecture, deserves an additional comment. 
Before the appearance of the Lorentz et al. paper jl^, the commonly accepted 
value of ph, based on an old calculation of Baym et al. ||l^, was some 50% 
higher. Consequently, the present results for /crust and A/dust are about two 
times smaller, than those obtained using the 'old' value of ph fl^ . 



5. SOLID CORES 

The appearance of crystal structure for p < ph is due to electromagnetic 
(Coulomb) interactions. Possible existence of solid cores in the liquid inte- 
rior of neutron stars would be due to specific features of the strong interactions 
between nucleons. At supranuclear density, localization of nucleons, character- 

, ^, |l^, |l|, 1^, H, results in 
the increase of the mean kinetic energy, Ackin- In order to make the crystal 
structure energetically preferred over the liquid one, the gain in the potential 
energy, Acpot, should more than compensate the increase of kinetic energy: 



istic of solid neutron star cores, discussed in |15 
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Acpot < ^Ackin- Calculations were based on variational approach: for an 
assumed nuclear hamiltonian, the structure of the ground state is usually de- 
termined from e{n\j, parameters) = minimum. Transition from liquid to solid 
takes place, if minimization on the crystalline (periodic) structure yields lower 
e than that obtained for homogeneous liquid. 

The possibility of solidification of pure neutron matter at sufficiently high 
density, due to strong repulsive short-range neutron-neutron interaction, was 
pointed out in the early 70-ties. However, subsequent calculations, based on 
more reliable solutions of the nuclear many-body problem, did not confirm these 
suggestions (see Q for the history of the neutron sohd). In the mid-seventies, 
several authors suggested existence of solid structures due to appearance of 
the neutral-pion-condensate ([^, @)- In all cases, solid structures were 
preferred because of the strong tensor component in the nucleon-nucleon in- 
teraction, which correlated nucleon spins with spatial directions. The specific 
crystalline structure depended on the assumed nucleon hamiltonian, while the 
spatial periodicities were simply related to the wave-length of the 7r° conden- 
sate, A^o. In the calculation, reported in the crystalline structure was a 
3-D cubic lattice of nucleons, with lattice spacing a = ^X^^o. In the paper of 
Takatsuka and Tamagaki dealing with pure neutron matter with tt" con- 
densate, the energy was minimized by a '1-D solid', consisting of alternating 
layers of neutrons, with opposite spin orientations, the distance between the 
layers being given by d = ^A^rO. It is more appropriate to call such a structure 
a 1-D liquid crystal (detailed review of the work of the Kyoto group on the 
crystal-like structures in pion-condensed neutron star matter can be found in 
20 ). In all cases, crystallization took place at p = 3 — Spo- 

A distinctly difercnt model of solid neutron star core was proposed by Kutschera 
and Wojcik |2^. As these authors show, for sufficiently low Xp < 0.05, pro- 
tons may be treated as 'impurities' in neutron matter, and their behavior is 
determined by their interaction with neutrons, Vnp. At the densities consid- 
ered, Vnp becomes more repulsive with increasing density, and this results in 
the decrease of neutron density around a p-impurity. For some nucleon-nucleon 
interactions, the crystalline structure appears at p > 3 — 4 po, with protons 
localized in the potential wells around the neutron density minima. 



6. DEFORMATION, ELASTIC STRAIN AND MOUNTAINS 

In contrast to the liquid phase of dense matter, the solid phase can support a 
non-vanishing shear and can be a site of elastic strain. We begin the present 
section with a brief summary of the basic concepts of the theory of elastic 
media. In what follows, for the sake of simplicity we will restrict ourselves to 
the standard, Newtonian version of the theory of elasticity . A consistent 
general-relativistic theory of elastic media can be formulated; this was done by 
Carter and Quintana p4[ |. 

The state of thermodynamic equilibrium (ground state at T = 0, if thermal 
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effects are neglected) corresponds to specific distribution of (equilibrium) posi- 
tions of ions within the neutron crust, r = (xi , X2 , x^). In what follows, we will 
neglect thermal effects. Such a state corresponds to the minimum of the energy 
density, e = eg = minimum, and to vanishing elastic strain. The deformation 
of an element of the crust with respect to the ground state configuration of ions 
implies a displacement of ions into new positions r' = r + u, where u = u(r) 
is the displacement vector. This is accompanied by the appearance of elastic 
strain (i.e., forces, which tend to return the system to the equilibrium state), 
and deformation energy density Ae = e — eg. Uniform translation, correspond- 
ing to r-independent displacement field u does not contribute to Ae, and the 
real (genuine) deformation is described by the (symmetric) strain tensor 

I ( du, duk\ , . 

'^^fe = o Ur- + 17- • (12) 



2 \ dxk dxi 

The non-zero elastic strain contributes to the stress tensor of matter, which for 
zero strain has the form —pSij , and which in the presence of elastic strain is 
given by 

= + crf^-" ■ (13) 

In the case of an isotropic solid, which should be a reasonable approximation 
to the macro scopically averaged neutron star crust (consisting presumably of 
small, disordered on macroscopic scale, anisotropic domains), the elastic prop- 
erties are determined by two elasticity parameters: the compression modulus 
and the shear modulus, ^. After the deformation, the volume of an ele- 
ment of matter changes according to dV' = (1 + Ujj)AV , where summation is 
understood over the repeated indices. The diagonal strain tensor, Uik = aSik, 
corresponds to pure compression. Compression term modifies the equation of 
state of the crust, adding an elastic contribution to the pressure p, Eq.(p^). 
The elastic shear contribution to the stress tensor reads 

a^]^'^^^- ^fi(2u,k-lujjS,k^ . (14) 

In the case of Ujj = we deal with a pure shear deformation of matter ele- 
ment, with crf^'^'^'' = 2/iUifc. In general, elastic strain involves both compression 
and shear. 

For a specific model of solid neutron star crust, the value of fi can be ex- 
pressed as a function of matter density and composition. For a bcc lattice 
structure below pnd one gets p5[ | 

2 

= 0.295Z^e^ni = 0.024 (^^^ ' pe , (15) 

where pe is the electron pressure. For p < pND we have p pe, so that ^ <C p. 
In view of the fact, that considered deformations will be small, the elastic effects 
can be always treated as small perturbation of a fluid (i.e. zero elastic strain) 



14 



P. Haensel 




9 10 11 12 13 



logio P (g cm 3) 



Fig. 4. — Shear modulus of neutron star crust versus density. Solid line - ground 
state of cold dense matter. Dotted line - accreted neutron star crust. 



configuration. For p > pnd, one has to take into account the presence of a 
neutron gas outside nuclei, but the inequality |crf^'^^'"| ^ p is still vaUd. The 
shear modulus of the neutron star crust, for the ground state composition of 
dense matter and for the accreted crust, is plotted versus density in Fig. 4. 

In the case of exotic nuclear shapes, predicted for the density interval ^Po < 
p < ph, numerical estimates of p, are not available, and in practice one is 
thus forced to use the extrapolation of the lower density (bcc coulomb lattice) 
results. 

In the case of solid pion-condensed core, shear modulus has been estimated 
in |l9|. In both cases one gets an estimate pcorc ~ 10~^ — 10~^ p, so that 
the elastic stresses of a hypothetical pion-condensed solid core could contribute 
significantly to the total stress tensor of dense matter. Actually, the perfect 
solid structure considered in Jl^ had a one-dimensional character (alternat- 
ing layers); however, the macroscopically averaged core should be expected to 
behave as an isotropic solid. In the case of bcc lattice of localized protons, con- 
sidered in p^ , we expect that pcoio << P, because of the small proton fraction 
in dense matter (the smallness of proton fraction was actually necessary for 
proton localization). 

The presence of elastic strain modifies the equation of mechanical equilibrium 
of neutron star. The equation of mechanical equilibrium of an element of 
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neutron star reads, for Newtonian gravitation, 

axj oxi 

where gt is the focal value of gravitational acceleration. The presence of af^^^^ 
allows for the permanent deformations of the neutron star crust. These defor- 
mations, supported by elastic stress, can lead to non- sphericity of non-rotating 
neutron star, and non-axiality of rotating neutron star. Such deformations - 
expected to be small (see below) - bear some resemblance to 'mountains' on 
surfaces of solid planets. The upper bound on the height of the 'mountains' on 
the neutron star crust will be discussed in the final part of this section. 

Another interesting possibility, stemming from the partially solid structure of 
neutron star, is that of precession. Precession is possible in the case of the non- 
equal principal moments of inertia, e.g., when I^x = lyy Izz- Precession of 
partially solid neutron stars, which produces a non-vanishing value of effective 
tri-axiality parameter, ecs, will be discussed in the subsequent section. 

Up to now, we restricted ourselves to Newtonian theory of elasticity. General 
relativistic theory of elasticity has been formulated in |Q. In its version suit- 
able for application to isotropic solids, relativistic theory involves two material 
parameters K and fx - which in the limiting case of weak gravity coincide with 
their Newtonian counterparts. The theory formulated in has been applied 
to fully relativistic description of partially solid, rotating neutron stars in . 

'Mountains' on the neutron star crust could be formed as a result of the 
neutron star crust-quakes. Such crust-quakes are expected to be due to cracks 
(disruption of the crystal structure) in the crust, implied by sufficiently large 
stresses (exceeding the critical values that the crust can support), and arising, 
e.g., during the slowing down of pulsar rotation. Let us consider a partially 
solid neutron star, rotating uniformly at the angular frequency Q, = 27r/P 
around the z - axis, which is assumed to be a principal axis of the moment 
of inertia tensor. Let us characterize the deformation of the neutron star by 
the asymmetry in the principal moments of inertia, {Ixx — lyy) 1 1 = where 
I = Izz- Such a rotating, non-axisymmetric neutron star, will be a source of 
pulsating gravitational radiation, of angular frequency ricR, = 251. The mean 
power of gravitational waves can be calculated using the formula derived in 
Section 16.6 of §, 

i^GR = f ^/W. (17) 

This formula can be used in order to confront the hypothesis of non- vanishing 
e with pulsar timing. Emission of gravitational radiation implies additional 
loss of the kinetic energy of the pulsar, and the slowing down of its rotation. 



according to ( i^kin = ^ -Bgr. In view of the fact, that i?kin = I^^, this 
V /GR . 

implies corresponding contribution to P, which will be denoted by (P)gr- 
Clearly, the slowing down, implied by gravitational radiation should be less 
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Table I. — Timing parameters and tmax for three selected radio pulsars 
Pulsar: 0437-47 1957+20 Crab (0531+21) 

P(nis) 5.76 1.61 33.4 

P (10-19 s s -1) 1.2 0.17 4.2 X 106 

d (kpc) 0.14 1.53 2.5 

emax 3 X 10-^ 2 X IQ-^ 1.2 x 10-^ 



that the observed one, [P] < Pobscrvod. This implies an observational 

V /GR 

upper bound on the value of e, 

^<^n.ax = 3xl0-(^)^(/^y ' 

where we assumed a 'canonical' value of / = 10"*^ g cm^. The values of e^ax for 
the Crab pulsar and for two selected millisecond pulsars, together with pulsar 
timing parameters, are given in Table ^. 

The choice of pulsars in Table | deserves a comment. The Crab pulsar was 
chosen as a good example of a young, seismologically active pulsar (glitches 
observed in the Crab pulsar timing are thought to result from neutron star 
crust-quakes). The remaining two pulsars are the millisecond ones; short rota- 
tion periods are preferred for an efficient radiation of gravitational waves, Eq. 
(0)- The PSR 0437-47 is very close, while PSR 1957+20 has the lowest value 
of P, and the second shortest period of all observed radio pulsars: its timing 
yields the most stringent upper bound on e. 



7. SHAPE AND ENERGY OF ROTATING, PARTIALLY SOLID 
NEUTRON STAR 

The existence of solid regions in the interior of neutron star (well established 
solid crust, speculative solid core) opens new possibilities for the structure 
and dynamics of rotating neutron stars. Partially solid neutron star can be 
axially non-symmetric, and can precess. In what follows, we will describe 
these properties of partially solid neutron stars with a simple 'one parameter' 
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model of Baym and Pines (hereafter referred to as BP). This model is 
Newtonian, and therefore enables one a very simple discussion of kinematics 
and energetics of rotating neutron star. Of course, neutron stars are relativistic 
stellar objects, with 2GM/Rc^ ^ 0.3, and the model of Baym and Pines is an 
approximation of reality. A consistent, general relativistic approach, based 
on the exact, relativistic theory of elasticity, was formulated by Carter and 
Quintana ||2^ (hereafter referred to as CQ). Connection between BP and CQ 
parameters will be discussed at the end of this section. 



7.1. Liquid neutron star 

Consider rigid rotation of a liquid neutron star around the z-axis, with angular 
velocity Q. Its stationary configuration will be axially symmetric with respect 
to z-axis which will, at the same time, be one of the three principal axes of 
the moment of inertia. Due to rotation, I^x — lyy < Izz- The deviation of 
the shape of the star with respect to the spherical one is described by the 
'oblateness parameter' 

e=^^, (19) 

where / = {Ixx + lyy + Izz)/-^- The total energy of the star is an even function 
of e, and for e ^ 1, and at fixed total angular momentum L is given by 

E = E, + ^ + Ae\ (20) 

^^zz 

where i?* is the total energy of a non-rotating configuration, and where terms 
higher than quadratic have been neglected. Stationary configuration corre- 
sponds to E = minimum^ at L = Izz^ = const. This yields the expression for 
the value of oblateness parameter of a liquid star, 

where /, is the moment of inertia of a non-rotating star. For the most rapid 
pulsars {P — 1.6 ms) the typical value of en is a few times 10^^ (assuming the 
mass AI ~ 1.4 M©) the exact value of eo depends both on the equation 
of state of dense matter, and on the neutron star mass. 



7.2. Rotation and elastic strain 



In general case, the solid regions of rotating neutron star are a site of elastic 
strain, and in particular, they are characterized by non- vanishing shear stresses. 
The elastic shear contributes to the local energy density of the solid region of 
neutron star, with Cshear = ^o'tkUik — fiUikUik- The contribution of shear 
stress to the total energy of the star is in the BP model represented by a 
simple expression, 

E'''^"^' = B{e - . (22) 
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The above formula represents a very important ingredient of the BP modeL The 
oblateness e corresponds to 'reference relaxed configuration' with zero shear 
(actually, with zero elastic strain) , consisting of the same number of baryons 
N. This configuration will be referred to as C. The most natural assumption 
is that C was the configuration of neutron star just after the solidification of 
the crust - so that its structure coincided, to a very good approximation, with 
that of a fluid star. The 'elastically relaxed' conflguration was rotating at a 
well defined Cl. By definition, E''^°''' = for n ^ il. 

The one-parameter formula of the BP model, which includes the effect of 
elastic shear on the rotating configuration, reads 

E ^ + — + Ae^ + B{e ^ ef , (23) 

where by definition of configuration C, 



dl 



AA \ d. 



(24) 

si=o 



Stationary, rotating configuration with a fixed L corresponds to E = minimum. 
This yields the equilibrium value of e. 

All effects of elasticity enter through the parameter B, another important pa- 
rameter is the oblateness of the 'relaxed reference configuration, e. Simple 
estimates of B for neutron star models with M 1.4 Mq with liquid interior 
yield B 10"^ A the effect of elasticity of the solid crust on the dynamics 
of rotating neutron star can thus be treated as a small perturbation. In view 
of i? ^ A it is suitable to introduce a small dimensionless 'rigidity' parameter, 
which measures the relative importance of elastic shear: b = B/A. According 
to ||2^, 6crust ~ 10~^. From Eq.([25|), the contribution of the elastic strain to 
the oblateness of rotating neutron star is given by (e)ciast — ^crust^. In the 
hypothetical case of a solid core in the pion-condensed neutron star interior, 
the effect is expected to be orders of magnitude larger, 6corc ^ 10^^ — 10^^ 
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General relativistic version of one-parameter formula for partially solid, ro- 
tating neutron star was derived by Carter and Quintana [^6| . Instead of making 
expansions in terms of the oblateness parameters e, e. Carter and Quintana use 
and Cl^ as expansion parameters. In the slow rotation limit, the Newtonian 
parameters A and B are related to the relativistic parameters P, and of 
ref.[|6|, by ^ 

^* " 2iA*+B) ■■ ^* " 2A{A + B) ' ^^^^ 
where /» is the moment of inertia of a spherical, non-rotating configuration. 
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8. OBLATENESS AND PRECESSION 



Even if originally fluid neutron star rotated around its symmetry axis, this 
situation might change after the star became partially solid. Let us consider 
rotating, partially solid neutron star. Let the instantaneous axis of rotation be 
defined by the unit vector n^. The deformation of neutron star results from 
the combined action of centrifugal forces and elastic stresses within the neutron 
star crust. Both effects can be treated as small perturbations. Let us switch- 
off deformation due to rotation; the resulting static configuration will still be 
deformed, due to elastic stresses within the neutron star crust. The symmetry 
axis of such a static but deformed 'reference' neutron star will be assumed to 
be given by the unit vector ng. In the case of partially solid neutron star the 
vectors no and ng, may differ: this might result from the previous action of a 
torque, which tilted the symmetry axis of the solid component of the rotating 
body with respect to the rotation axis. Such torques, which could have acted 
on neutron star crust in the past, could be exerted by the neutron star magnetic 
field, by the tidal forces from a close companion, or by a passing-by neutron 
star. Tilting could also result from a series of 'asymmetric' crust-quakes. In 
the general case, the unit vectors n^ and Uq will not be fixed in space, but 
will rotate around the fixed direction determined by the angular momentum 
vector L. The rotation of n^ around L corresponds to well known precession 
of a solid body, typical for the case when the rotation axis does not coincide 
with any of the three principal axes of the body moment of inertia. Keeping 
only terms linear in oblateness parameters, we may write the inertia tensor of 
rotating neutron star as pR M 



(1 



1 



1 



be)S, 



(27) 



The angular momentum of the star is given by Li — lij^j. The three vectors 
L, no and no are coplanar pl[ . 

For a completely liquid star, we would have 6 = 0, and we could take n^ as 
z-axis. This would give 



^liquid J 



1 3 



(28) 



which would imply 



rliquid rliquid 

XX yy 

^liquid 



(29) 



In the above equations, is the moment of inertia of a spherical, non-rotating 
neutron star. 

In the linear approximation, the deformations described by the oblateness 
parameters en and be can be treated independently. From the Euler equations. 
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one finds that the angular velocity vector (as calculated in the solid body ref- 
erence system) precesses around the direction of no, with precession frequency 

Oprec - . (30) 

The free precessional frequency of a partially solid neutron star is thus deter- 
mined by internal stellar structure and distribution of internal elastic strain. 

While it might appear, that the free precession of partially solid neutron 
stars is a rather common phenomenon, observations of radio pulsars give only 
one (and debatable) example of the possible presence of precession in the pulsar 
timing. Analysing timing data for the Crab pulsar, spanning the time interval 
of five years, Lyne et al. ||2^ find 'quasi-sinusoidal' deviations from the regular 
pulsar slow-down, with an approximate period of approximately 20 months. It 
should be stressed, that originally Lyne et al. preferred the interpretation of 
it in terms of the so called 'Tkachenko oscillations' of the superfluid interior 
of neutron star. Further quantitative considerations indicated, however, the 
relevance of the solid interior of neutron star for explaining the specific period 
of modulations IpOf . Precession of the solid crust remains thus one of possible 
explanations |3C[| . Applying our precession model for explaining observed quasi- 
sinusoidal deviations, we would get (6e)crab ~ 10^^". 



9. GRAVITATIONAL RADIATION FROM PRECESSING PUL- 
SARS 



9.1. Mean power of gravitational radiation 

Precession of a deformed pulsar implies radiation of gravitational waves. The 
mean power of gravitational radiation is calculated using the general formula, 

where is traceless mass quadrupole moment in laboratory (observer's) sys- 
tem [0. In the lowest order in e's, terms proportional to en do not contribute 
to the relevant time derivatives, and therefore when calculating Xij we may put 
en = 0. The calculation gives 0]: 

Eon = ~ih~h)l^o^'0' , (32) 

where 9 is the 'wobble angle' (the angle between the reference ('sohd body') 
symmetry axis no and (fixed in laboratory, inertial reference frame) direction 
of L; we assumed that 6* -C 1. In Eq. (|3^), Ii and I2 are principal moments 
of inertia in the neutron star (body) reference system (/i = I^x, I3 = Izz)- 
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The frequency of the gravitational waves is approximately equal to the pulsar 
frequency, fl, namely ^on = + ^be^ In our case, we have 

h-h^ ^beh , (33) 
so that the formula for -Eqr can be written in the form 

Eon - f^l'^'^s , (34) 
where the effective tri-axiality is 

ecff = Ibee , (35) 

o 

and I = I^,. 



9.2. Estimating b, e, and 9 

For a neutron star of mass M 1.4 Mq the existing estimates give fcciust ^ 10~^ 
. In the hypothetical case of a more massive neutron star with a solid pion 
condensed core, a very uncertain estimate is 5corc ~ 10^'^ — 10^^ [2^ . 

Let us first calculate the value of e as a function of angular frequency of radio 
pulsar. Using Eq.(|2l]), we get 



'-'^'-^^ ' (lO^gcmO (l053erg) (lOms) ' ^^^^ 

Neutron star, born with some initial period of rotation, Pi„it , slows down during 
its life as a radio pulsar. In Table || we give numerical estimates of initial 
and present oblateness parameters, einit and e for the most rapid millisecond 
pulsar, PSR 1937+21, and for the Crab pulsar. In both cases we assumed 
/ — W^^ g cm^ and A = 10^^ erg. Initial period of the Crab pulsar was 
calculated using the magnetic dipole formula for pulsar slowdown In the 
case of the PSR 1937+21, we considered two cases. The milhsecond pulsars 
are believed to be previously radio-dead, old pulsars (P ^ seconds) , revived 
by accretion of matter in close binary system. In the first case, we assumed 
that the old pulsar (predecessor of PSR 1937+21) had in the beginning of the 
accretion era small deformation, equal to eo at the beginning of this era, when 
P ~ Pinit = 1 s. We assumed additionally, that this initial configuration was 
elastically relaxed (no elastic strain) . Notice, that in the case of the Crab pulsar 
the initial configuration was hot and liquid, while in the case of PSR 1937-21 
the pulsar was assumed to be sufficiently old to be elastically relaxed. 

In principle, in the case of PSR 1937+21, we could contemplate a rather 
extremal situation, when the pulsar keeps the initial oblateness from its 'first 
birth'. Such a situation, which would require a very high resilience of the 
stellar solid, would lead to a Crab-like initial oblateness, which is indicated 
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Initial and present oblateness parameters of two pulsars 



Pulsar: 1937+21 Crab (0531+21) 



P 1.56 ms 33 ms 

Pinit 1 s (10 ms) 17 ms 

einit 10-7 (10-4) 3 X 10-4 

en 4 X 10-2 9 x 10-^ 



within the brackets. Notice however, that during the accretion era the old 
pulsar had to accrete some '^0.1 Mq, in order to spin-up to P ~ 1 ms. In view 
of this, the 'old crust', with large initial oblateness, was presumably 'molten' 
by compression, and transformed into the outer layer of the liquid interior. It 
seems thus, that in order to keep the 'first birth' oblateness, the millisecond 
pulsar would have to contain a solid core. 

The value of the 'wobble angle', 0, is limited by the critical strain, which 
the solid crust can support. The crust breaks down (fractures) if the shear 
exceeds some critical value. In terms of oblateness parameters, the condition 
for disruption of the solid is |e — e| > Ocrit, where Gdit is called 'critical angle' 
(notice that shear deformation changes angles within the crystalline structures) 
[32| . A non-zero wobble angle implies an elastic strain within the solid region 
of the star. For typical terrestrial solids, Gcrit ~ 10-^ — IQ-*. This would 
imply e ~ CQ (a stringent condition in the case of millisecond pulsars). Such 
a low Gcrit would limit the wobble angle to ^^<10-"^ [Q. Let us notice, that 
for a perfect Coulomb lattice - which seems to be better suited as a model 
of neutron star crust - calculations give 8crit ^ 10~^ — 10~^, which would 
imply e ~ Cinit, and no limitations on 6 ! (An interesting discussion of 'elastic 
strength' of neutron star's crust, referring also to terrestrial experiments on 
'near Coulomb' metallic microcrystals, can be found in |]3^ .) However, neutron 
star crust should expected to be very far from the 'perfect crystal' structure, 
and therefore we may expect a significantly lower value of Gcrit- In any case, 
it is difficult to point out a mechanism, which could build up a wobble angle 

0<i m. 



Table III. — 
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Pulsar: 1937+21 1957+20 Crab(0531+21) 



P(ms) 


1.56 




1.61 




33 


P (10-19 s s-i) 


1.0 




0.17 




4.2 X 10^ 


^max 


6 X IQ- 


-9 


2 X 10^ 


-9 


1.2 X 10-3 


en 


4 X IQ- 


-2 


4 X IQ- 


-2 


9 X 10-5 


eeff - case a 


1 X 10- 


10 


1 X 10- 


10 


3 X 10-" 


Eoff - case b 


4 X 10" 


13 


4 X 10" 


13 


3 X 10-12 



10. GRAVITATIONAL RADIATION FROM PRECESSING PUL- 
SARS AND PULSAR TIMING 

The slowing-down of the pulsar rotation due to gravitational radiation is limited 
by the measured value of P, so that (P)gr < Pobs- This yields observational 
constraint on the effective tri-axiality of pulsar, 



eeff < emax = 3 X 10 ^ 



P \ ^ P 



1ms/ \ 10-19 



eeff = Ibee. (37) 

o 

We will consider two limiting cases, of a small and large critical strain, respec- 
tively. Numerical estimates for these two limiting cases, based on the timing 



data for two millisecond pulsars and the Crab pulsar, are given in Table [II. 
In all cases, we assumed rigidity parameter b = 10"^. For the scenario corre- 
sponding to 'case a', we use a conservative assumption of = 10-'^. In view 
of the assumed small value of critical strain, we have in this cas e ~ eo. Due 
to rapid present rotation of millisecond pulsars, their effective tri-axialities are, 
in 'case a', some three orders of magnitude higher than this of the Crab pul- 
sar. The situation changes in 'case b' (large critical strain). Then, we could 
contemplate the value of 9 for millisecond pulsars, which could be as large as 
1; however, due to the large critical strain, we should expect e = einit, corre- 
sponding to Pinit = 1 s, which reduces the effective tri-axiality by some three 
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orders of magnitude as compared to the 'case a'. For the Crab pulsar we put 
e = Cinit) but in view of its short age (941 years) we made a conservative choice 
oi = 10~^ (no sufficient time to built a large tilt of the 'solid body axis' with 
respect to rotation axis). Under these assumptions, we get in 'case b' similar 
values of effective tri-axiality for the most rapid millisecond pulsars and the 
Crab pulsar. 

The presence of a sufficiently large solid core would imply some difficulties in 
confrontation with pulsar timing. In such a (highly speculative) case, we would 
expect b ~ bcore ~ 10~^ — 10~^. Even with a rather small = 10~^, the large 
value of rigidity parameter would lead to effective tri-axiality largely exceed- 
ing the observational upper bound for millisecond pulsars. Moreover, as the 
critical strain for the solid core may be expected to be very large, the possible 
'wobble angle' that could be built during star evolution could, in principle, be 
large. Confronting these theoretical speculations with pulsar timing leads to 
conclusion, that if millisecond pulsars contain a solid core, then their wobble 
angles are extremely small. A more restrictive statement is that millisecond 
pulsars are just unlikely to contain a massive solid core. 



11. DETECTABILITY 



The radio pulsars with a non-zero tri-axiality could be promising sources of 
continuous gravitational radiation. The most attractive in this respect are 
the nearby millisecond pulsars with eoff ~ 10~^ [Q. Consider two selected 
millisecond pulsars, PSR 1957-^20 and PSR 0437-47. The advantage of the 
first one is its very short period, P = 1.6 ms (the distance being d = 1530 pc). 
On the other hand, the second millisecond pulsar is much closer, at d = 140 pc. 

Millisecond pulsars with non-zero tri-axiality would be interesting sources 
of continuous, nearly monochromatic gravitational radiation, provided some 
additional conditions are satisfied. As we could see, the fact that millisecond 
pulsars are such precise clocks (P < 10^^^ s s"^ !), imposes rather stringent 
constraints on the maximum value of Coff. The timing of the millisecond pul- 
sars tells us, that the possible amplitude of the gravitational waves should be 
sufficiently small to be consistent with observations. As a consequence, the 
detectability of gravitational waves from the nearby millisecond pulsars by the 
forthcoming generation of the LIGO and VIRGO detectors would be possible 
only if the periodic signal could be integrated over the time of a few years Q . 
This would require a sufficiently long lifetime of the precession or wobble. 

In the approximation, in which the coupling of the solid and liquid regions 
of neutron star are neglected, liquid interior of neutron star cannot partici- 
pate in stellar precession. However, such a coupling exists, and detectability of 
continuous gravitational radiation requires that the damping times connected 
with the liquid-solid coupling be substantially longer than a few years ||3^ . 
Finally, the 'wobble' is damped by the gravitational radiation itself (see Exer- 
cise 16.13 in iQ]). Explicit calculation shows that the characteristic timescale 
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^wobble ~ lQ^{he/lQ-^)-^{Pll ms)-'^ years (see, e.g., ||): the detectability 
requires, that t^o^j^i;, be substantiaUy longer than a few years. Our discussion 
of the damping timescales indicates thus, that the situation in which a non- 
zero tri-axiahty is due to sufficiently large, genuine 'mountains' on the neutron 
star crust, is the most attractive one. This conclusion seems to be in a sur- 
prising harmony with the location of Les Houches School, dominated by the 
magnificent massif du Mont Blanc, so admired by all of us during the coffee 
breaks. 
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